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We outline a method to infer the global history of star formation in galaxies with 
input only from absorption-line observations of quasars. The application of the 
method to existing data leads to the conclusion that most stars formed at rela- 
tively low redshifts (z < 2). We combine the global rate of star formation with 
stellar population synthesis models to compute the mean comoving emissivity and 
mean intensity of background radiation from far-UV to far-IR wavelengths. These 
predictions are consistent with all the available measurements and observational 
limits, including recent results from HST and COBE. 

1 Overview 

This article concerns the evolution of, and relations between, various large- 
scale average properties of the population of galaxies as a whole. It is often 
convenient to express these "global" properties as mean comoving densities 
and to normalize them to the present closure density. We are particularly 
interested in the comoving densities of stars, gas, metals, and dust within 
galaxies, which we denote respectively by fl s , Sl g , £l m , and fid- The last three 
of these are meant to refer to the interstellar media (ISM) of galaxies, exclusive 
of the intergalactic medium (IGM), although in practice such a distinction 
may only be approximate. As defined here, O m includes metals in both the 
gas and solid (i.e., dust) phases of the ISM. It is usually more informative to 
reexpress Q m and fid in terms of the mean metallicity and mean dust-to-gas 
ratio, Z = fl m /fl g and D/G = fld/fl g - It is clear that all of these properties 
are related in the sense that, as new stars form, fl s will increase, while, in most 
cases, fl g will decrease, and Z and D/G will increase. One of our goals is to 
quantify such relations through the equations of "cosmic chemical evolution" . 

Until recently, there were no emission-based estimates of the global rate of 
star formation Q s at z > 0.3. The reason for this is that samples of galaxies se- 
lected by emission become progressively incomplete and include only brighter 
objects at higher redshifts. In contrast, samples of galaxies selected by absorp- 
tion against background quasars do not suffer from this bias. Such observations 
are exquisitely sensitive to small column densities of absorbing or scattering 
particles. In principle at least, they enable us to estimate fl g , fi m , and even 
fid as functions of rcdshift. From these and the equations of cosmic chemical 



evolution, we can then infer the global rate of star formation fl s . It is amusing 
to note that this idealistic program does not require the detection of a single 
stellar photon! Furthermore, if we are confident (or foolish) enough, we can 
combine our estimates of fl s with stellar population synthesis models to com- 
pute the mean comoving emissivity £ v and the mean intensity of background 
radiation 3 V . One might then claim to have predicted the "emission history" 
of the universe from its "absorption history" . This article describes a first at- 
tempt by Yichuan Pei, Stephane Chariot, and the author to carry out such a 
program; a complete account of our work is given in references 1 and 2. Some 
related material can be found in references 3-6. 

2 Absorption-Line Systems 

Before proceeding, it is worth recalling some facts about the statistics of 
absorption-line systems. Let f(N x ,z) be the column density distribution of 
particles of any type x that absorb or scatter light. These might, for example, 
be hydrogen atoms (x = HI), metal ions [x = m), or dust grains (x = d). By 
definition, Ho(l + z) 3 \dt/dz\f(N x , z)dN x dz is the mean number of absorption- 
line systems with column densities of x between N x and N x + dN x and redshifts 
between z and z + dz along the lines of sight to randomly selected background 
quasars. These lines of sight are very narrow (much less than a parsec across) 
and pierce the absorption-line systems at random angles and impact parame- 
ters. One can show that the mean comoving density of x is given by 

^) = ^r L dN x f(N x ,z)N x , (1) 

where m x is the mass of a single particle (atom, ion, or grain). Equation (1) 
plays a central role in this subject. It enables us to estimate the mean comoving 
densities of many quantities of interest without knowing anything about the 
structure of the absorption-line systems. In particular, we do not need to 
know their sizes or shapes, whether they are smooth or clumpy, and so forth. 
A corollary of equation (1) is that the global metallicity, Z = O m /O s , is 
given simply by an average over the metallicities of individual absorption-line 
systems weighted by their gas column densities. 

The absorption-line systems of most interest in the present context are the 
damped Lya (DLA) systems. It is widely believed that they trace the ISM of 
galaxies and protogalaxies and are the principal sites of star formation in the 
universe. There are excellent reasons to adopt this as a working hypothesis. 
First, the DLA systems have, by definition, iVni > 10 20 fp -2 , which is just 
below an apparent threshold for the onset of star formationQ. Second, the DLA 
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systems-contain at least 80% of the HI in the universe and appear to be mostly 
neutrala. The other absorption-line systems, those with Nm < 10 20 cm~ 2 , 
probably contain more gas in total than the DLA systems, but this must be 
diffuse and mostly ionized. In the following, we regard non-DLA systems as 
belonging to the IGM, even though some of them might actually be located 
in the outer halos of galaxies. This distinction - between the mostly-neutral 
ISM, where stars form, and the mostly-ionized IGM, where they do not - is 
clearly valid at the present epoch. Thus, the DLA systems are often referred 
to as DLA galaxies. It will be interesting to see exactly which types of galaxies 
they represent, but as we have already emphasized, this issue does not affect 
any of the global properties derived from equation (1). 

The sample of known DLA galaxies now includes about 80 objects!. They 
are distributed over a wide range in redshift, < z < 4, although, as a con- 
sequence of selection effects, most of them are confined to the narrower range 
2 z < 3. From observations of DLA galaxies in various subsets of this sample 
and comparisons with present-day galaxies, the following trends have emerged. 
The mean comoving density of HI decreases by almost an order of magnitude, 
from H i » (1 - 2) x lO" 3 ^" 1 w 3 to fl m w 2 x lO^h' 1 at z = [with 
h = Hq/(100 km s _1 Mpc -1 )] Hfla. It is possible that fijji increases between 
z w 4 and z w 3, but the evidence for this is weakla. The mean metallicity 
increases by about an order pfpiaaagnitude, from Z w 0.1 Zq or slightly less at 
z ~ 2 to Z w Zq at z = EBtil. The mean dust-to-gas ratio increases by a 
similar factor, while the mean dust-to-metala iratip remains roughly constant 
at about the present value in the local ISM t3'lliHl3. Thesej^suJis are entirely 
consistent with the recent Keck observations by Lu et al The,-abun- 
dances of H2 and CO appear to be much lower at z > 2 than at z — OtB. As 
a consequence of the relatively small samples involved, most of the numbers 
quoted here are uncertain by factors of 1.5 or more. 



3 Cosmic Chemical Evolution 

The global properties defined above are governed by a set of coupled equations, 
which are sometimes referred to as the equations of cosmic chemical evolution. 
In the approximation of instantaneous recycling (and Z « 1), they take the 
form 

j t (n g + n s ) = n f , (2) 

^(zn g ) + (z-y)jn s = z f n f , (3) 
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where y is the IMF-averaged yield. Equations (2) and (3) are strictly valid only 
when all galaxies evolve in the same way; otherwise, they should be regarded as 
approximations. The "source" terms on the right-hand sides of the equations 
allow for the exchange of material between the ISM of galaxies and the IGM; 
they represent the inflow or outflow of gas with metallicity Zf at a rate f2/. 
To illustrate a range of possibilities, we consider three types of evolution: a 
closed-box model (fi/ = 0), a model with inflow of metal- free gas (0/ = +v(l a , 
Zf — 0), and a model with outflow of metal-enriched gas ((if = —v(l s , Zf = 
Z) . Our inflow and outflow models are direct analogs of the standard modelsjQf 
chemical evolution in the disk and spheroid components of the Milky WayEjO. 
We fix the yield y in each model by requiring Z = Zq at z = 0. Then the 
only adjustable parameters are the "initial" comoving density of gas in galaxies 
Q goo (in practice, the value of Q g at z > 4) and the relative inflow or outflow 
rate v. 

To complete the specification of the models, we make two other approxima- 
tions, both motivated by the observations summarized in the previous section. 
(1) We neglect any ionized or molecular gas in the ISM of galaxies and set 
fig = 1.3fiHi (to account for He). (2) We assume that just over half of the 
metals in the ISM are depleted onto dust grains and set D/G = O.QZ. The 
models are designed to reproduce (as input) the observed decrease in the mean 
comoving density of HI between z « 3 and z = 0. The only subtlety here is 
that the observed values of Ohi tend to underestimate the true values as_a 
consequence of the obscuration of quasars by dust in foreground galaxies 113. 
We make a self-consistent correction for this bias in the models by linking 
the obscuration of quasars to the chemical enrichment of galaxies. It is worth 
noting that, while this correction has a substantial effect on flm, especially at 
z ~ 1, it does not entail large numbers of "missing" quasars (only ^20% at 
z = 2 and ~40% at z = 4). The models reproduce (as output) the observed 
increase in the mean metallicity between z w 2 and z — without any fine 
tuning of the parameters £l goo and v. The reason for this is that most of the 
star formation and hence most of the metal production occur at z < 2. 

Figure I shows the evolution of the comoving rate of metal production 
p z in the models. This is given by p z = yip, with tp = (1 — i?) _1 p s and 
p s = (3i?Q /8ttG)(I s , where R ~ 0.3 is the returned fraction. The predicted 
rates have maxima at I < z < 2 and decline rapidly at lower redshifts. Fig- 
ure I also shows estimates of, and lo wer linn its on, p z from recent ground-based 
surveys and the Hubble Deep Field Oocj'Ej. These are proxies for global Ha 
and UV emissivities based on the closepCorrespondence between UV emission 
and metal production in massive stars Uc3. Evidently, the predicted and ob- 
served rates are in broad qualitative, and even some quantitative, agreement 
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Figure 1: Comoving rate of metal production p z as a function of redshift z (for h = 0.5, 
qo = 0.5, and A = 0). The curves are from the closed-box (C), inflow (I), and outflow 
(O) models with f2 goo = 4 X 10 — 3 h~ 1 and v = 0.5 (see Figure 1 of reference 1). The data 
points and lower limits represent global Ho and UV emissivities from ground-based and HST 
surveys (see Figure 9 of reference 22). 



(given the uncertainties in both). This is remarkable because the models were 
constructed only with absorption-line systems in mind, not the emissivities rep- 
resented in Figure 1. We have also combined our chemical evolution models 
with stellar population synthesis models to compute directly the mean comov- 
ing emissivity £„ at wavelengths from 10 _1 /im to 10 3 /im and, by an integration 
over redshift, the corresponding mean intensity of background radiation J„. 
These calculations include a self-consistent treatment of the absorption and 
reradiation of starlight by the dust within galaxies. The same models shown 
in Figure 1 also predict a far-IR/sub-mm background in nice agreement with 
a tentative detection based on COBE dataE3. 
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